EMPG-02-23 
HWM-02-38 
hep-th/02l20T6 



On the existence of minima in the Skyrme model 



B. J. Schroers 
Department of Mathematics, Heriot-Watt University, 
Riccarton, Edinburgh EH14 4 AS, UK 
e-mail berndOma .hw.ac.uk 

20 November 2002 
Abstract 

Well-separated Skyrme solitons of arbitrary degree attract after a suitable relative 
rotation in space and iso-space, provided the orders of the solitons' leading multipoles 
do not differ by more than two. I summarise the derivation of this result, obtained 
jointly with Manton and Singer, and discuss to what extent its combination with 
earlier results of Esteban allows one to deduce the existence of minima of the Skyrme 
energy functional. 
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1 Introduction 



The Skyrme model is a non-linear field theory in which nuclei are modelled by quantum 
states of classical topological solitons . In this talk I will use the term topological soliton 
for minimal energy solutions in a field theory which have an associated integer degree or 
topological charge. The degree is conserved for topological reasons and ensures the stability 
of the soliton. In the Skyrme model static configurations are maps 

U : M 3 -> SU(2). (1.1) 

Points in M 3 will be denoted by x = {x\,X2,x%) and I write r for the Euclidean length 
|x| = \/ xf + x\ + x\. Physically, U combines the pion fields 7ri, 7^,773 and the a-field 

U(x) = a(x) +iir a (x)r a , (1.2) 
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where r a , a = 1, 2, 3, are the Pauli matrices and the constraint cr 2 + irf + 7r| + vr| = 1 is 
imposed. Here and in the following, summation over repeated indices is always implied. 
The Skyrme energy functional is 

E[U)=-J d 3 x Qtr(L J L 4 ) + ^tr([L„L J ][L„L J ])^ , (1.3) 

where 

Li = U^d t U (1.4) 

and di = d/dxi with % — 1, 2, 3. The Euler-Lagrange equation can be expressed in terms of 
the modified currents 

Li = Li — ^[Lj, [Lj, Li}} (1.5) 

and reads 

diU = 0. (1.6) 

It was already noted by Skyrme that the finite-energy requirement means that the Skyrme 
fields have to tend to a constant at spatial infinity, thus becoming effectively maps from S 3 
to SU (2) ~ S 3 with an associated integer degree. Skyrme identified the degree physically 
with the baryon number. The first rigorous proof that for finite-energy Skyrme configura- 
tion the degree 

deg[U) = -— J d 3 xe ijk tr(L i L j L k ) (1.7) 

is an integer was only given fairly recently in 0. As a result the configuration space 

C = {U : M 3 -> SU(2) | E[U] < oo} (1.8) 

is partitioned into sectors Ck consisting of all finite energy configurations of degree k. 
Faddeev's bound [Q 

E[U] > 12n 2 \k\ (1.9) 

(the strict inequality follows from ||]) implies the existence of infima 

I k = M{E[U] | U g C k }. (1.10) 

The central question of this talk is whether the infima are attained, i.e. whether minima 
exist in all sectors of the Skyrme model. In the following I call minimal energy solutions 
of non-vanishing degree Skyrme solitons. 

To illustrate quite how little is known rigorously about the existence of Skyrme solitons 
note that the hedgehog ansatz 

U H (x) = exp(if(r)x a r a ) (1.11) 
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with the boundary condition /(0) = 7r and /(oo) = leads to a minimisation problem for 
the profile function /, which is known to have a unique solution ||, called the Skyrmion. It 
is also known |J that the Skyrme energy functional has a minimum in the sector C\. How- 
ever, it is still not known whether the minimum is of the spherically symmetric hedgehog 
form ( |mi) . 

Physically, we may think of Skyrme solitons of degree k > 1 as bound states of Skyrmions. 
One therefore expects the questions of whether such bound states exist to be related to the 
existence of attractive forces in the Skyrme model. This expectation is borne out by the 
analytical work of Esteban in ||. Esteban showed that for a suitable class functions 

h<h + h-i (1.12) 

for all k, I G Z. Physically this result is only sufficient to ensure threshold bound states - it 
would be satisfied if for example Ik = C\k\ for some constant C. Esteban went on to show 
that a minimum exist in provided one assumes the strict inequality 

J fc </, + J fc _, (1.13) 

for all k e Z - {0, ±1} and I e Z - {0, k} in the range \l\ + \k - l\ < y/2\k\. This result 
will be referred to as Esteban's theorem in the following discussion. In it was shown 
that the result still holds if one widens the class of allowed functions, but the inequality 
( 1.13|) remains a necessary assumption in the proof. In the cases where the infima l\ and 



Ik-i are attained by Skyrme solitons, the inequality (|1.13|) is equivalent to the existence of 



attractive forces between those Skyrme solitons. In [Q] the existence of attractive forces is 
studied in detail for Skyrme solitons obeying certain regularity assumptions. In the second 
part of this talk I review the arguments and results of that paper and in the third part I 
discuss their relation to Esteban's result. 

Before plunging into the technical analysis I should point out that much is known numer- 
ically about Skyrme solitons. Numerical searches aided by analytical ansatze suggest the 
existence of Skyrme solitons for 1 < k < 22 || |9|, [HJ. However, it is also worth stressing 



that the existence of classical minima does not ensure the existence of a quantum bound 
state. This is relevant in the physical application of the Skyrme model because nuclei do 
not exist for arbitrarily large baryon number. The study of the deuteron as a quantum state 



of the k = 2 toroidal Skyrme soliton in [[LI]] furthermore shows that when such a quantum 
bound state exists the quantum mechanical matter distribution may be quite different from 
that of the underlying classical soliton. 



2 The interaction energy of Skyrme solitons 
2.1 Asymptotics of a Skyrme soliton 

The goal of this subsection is to study the behaviour of Skyrme solitons near spatial infinity. 
The key step is to think of the 2-sphere at infinity as a boundary of M 3 and to show that 
the Skyrme equation is regular there. It then follows from a unique continuation argument 
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that non-trivial solutions of the Skyrme equation have a non-trivial expansion near infinity. 
In other words, Skyrme solitons necessarily have a non-trivial large r expansion in powers 
of 1/r, possibly combined with lnr. 

To explain the basic ideas, consider a finite point Xq and assume (if necessary redefining 
U(x) — > U(xq)~ x U(x) ) that U(xq) = 1. Then we have the expansion 

U(x) = l + u{x) (2.1) 

in a neighbourhood of Xq, with the 2x2 complex matrix u satisfying the algebraic con- 
straints 

u + vt + uv! = 0, tr(u) + det(u) = 0. (2.2) 

After lengthy algebra, the Skyrme equation can be rewritten 

P(u, du, d 2 u) = Q(u, du) + F(u, du), (2.3) 

where P is linear in the second derivatives d 2 u and contains the Laplace operator Am. The 
functions Q and F crucially do not depend on d 2 u and are, respectively, quadratic and of 
degree four in u. The point of writing the Skyrme equation in this way is that the operator 

f^P(u,du,d 2 f) (2.4) 

is linear and elliptic near xq. If we assume that both u and the currents Lj are Holder 
continuous, we deduce by elliptic regularity that a solution of fl2.3p is twice differentiate. 
Bootstrapping further, one concludes that any solution is in fact smooth (i.e. C°°). The 
key idea in this procedure is to re-interpret the non-linear Skyrme equation as an elliptic 
linear equation with coefficients depending on u and du. 

In order to apply the same idea at spatial infinity, we treat the 2-sphere at spatial infinity 
as a boundary of M 3 and introduce coordinates (s,0,ip), where s = 1/r and 9 and tp are 
the usual spherical coordinates. Then defining 

A = rdi = ~di (2.5) 

s 

the Euclidean Laplacian takes the form 

A = s 2 A b (2.6) 

with 

A, = s 2 d 2 + A„ (2.7) 

where A w is the Laplacian of the unit 2-sphere. Now we write U = 1 + u for large r (that 
is, for small positive s) and obtain a '^'-version of the Skyrme equation 

P b (u, Du, D 2 u) = Q b (u, Du) + s 2 F b (u, Du). (2.8) 
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This equation can now be analysed using tools from the theory of so-called 6-differential 



operators, see |L2], pL3| . In particular one can use the boundary conditions that U tends to 
one at oo and that the currents Lj vanish there to solve fl2.8|) iteratively, thus producing 
an expansion of the solution in powers of s — 1/r, with possible factors lnr in higher order 
terms. At this stage we have an asymptotic expansion of Skyrme solitons, but do not yet 
know whether there is a leading non-zero term in the expansion. This can be established 



with the help of a unique continuation theorem |TJJ], suitably adapted for our purposes 



The upshot of this chain of arguments, described in detail in jTJ, is that every Skyrme 
soliton has a leading Lie-algebra valued multipole field 

M 

2M + l^ Mm r 1 



u M {x)=ir a ^ — — -Q Mm , (2.9) 



m=-M 

where Yum are the usual spherical harmonics on S 2 . In accordance with the usual nomen- 
clature we refer to the multipole field (|2.9|) as a 2 M -pole and call M the order of the 
multipole. The leading multipole moments Q1i m are independent of the location of the 
Skyrme soliton, and are acted on naturally by rotations and iso-rotations. They are crucial 
for the calculations in the following section. Note that it was already shown in [17] that 
the leading multipole cannot be a monopole. The Skyrmion field ( |1.11|) is known to have 
an iso-triplet of dipoles as leading multipoles. For other Skyrme solitons, too, the leading 
multipoles have been investigated to some extent. The leading multipole of the largest 
known order is an octupole (i.e. multipole order 3) which is believed to arise in the charge 
seven icosahedral Skyrme soliton || 

2.2 Interaction energy of two scalar multipoles 

As a technical preparation for the computation of the interaction energy of two Skyrme 
solitons we derive a formula for the interaction energy of two scalar multipoles. In our 



description of multipoles we use the conventions of [D| throughout. Consider the field of 
a 2 M -pole centred at X + = (0, 0,i?/2), where R > 0. In terms of spherical coordinates 
(9 + ,<f + ) centred at X + it reads 

MX) - (2M + 1) ^ QMm \ x -XA^ (2 ' 10) 



m=-M 



Similarly define X_ = (0, 0, —R/2) and let (6L, y9_) be spherical coordinates centred at X. 
A 2 w -pole field centred at X- = (0, 0, -R/2) has the form 



N 



9n{x) 



4tt ^ ~ Y Nn {0-,yJ) 

Y) 2^ QNn \ x -X\ N + 1 ' ^ ' 



{2N + _, 

v ' n =-N 



A natural measure for the interaction energy between harmonic function in the upper and 
those in the lower half plane is 

V[f,g] = [ d Xl dx 2 (gd 3 f - fd 3 g). (2.12) 
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The computation of this interaction energy for the multipole fields fu and is not easy 
because the combined field of the two multipoles only has cylindrical symmetry about the 
a;3-axis and not the full rotational symmetry of each of the multipoles. As explained in ||7J, 
Fourier transformation in the Xia^-plane turns out to be an efficient tool for evaluating the 
integral in ( j2.12j) . Assuming without loss of generality that M < N, the answer is 

_ (4tt) 2 (M + AQK-1)^ 

V V M 'M-R»+™ v / (2M + l)(2iV + l) = X 

m ^ M y/(M - m)\(M + m)\^(N - m)\(N + m)\ ' 



For us, the most important feature of the formula (|2.13|) is that, for non- vanishing multi- 
pole moments, the interaction energy can always be made non-zero by a suitable rotation 
of the multipole of the highest order, in our case N (if M = N it does not matter which of 
the multipoles gets rotated). By definition the multipole components Qnh are the compo- 
nents of an element Q of the (2N + l)-dimensional irreducible unitary representation Wn 
of SO (3), with G G SO (3) acting via Qnu | — ► J2n'=-N Un n >(G)QNn'- Think of the pairing 
of the multipole moments in (|2.13|) linear form 

F Q :W N ^R, Q^ Y QMmQNm 

w y/(M - mjUM + m)\y/(N - m)\(N + m)\ 

Then it follows from the irreducibility of Wn that U N (G)Q cannot lie in the kernel of Fq 
for all G. Therefore Fq(U n (G)Q) ^ for some G. Hence the interaction energy ( |2.13| ) is 
non-vanishing after rotating Q with that G. 



2.3 Interaction energy of well-separated Skyrme solitons 

Consider now two Skyrme solitons and of degrees k and I. Since the total energies 
of both and are finite there must be balls B\ and B 2 in M. 3 so that most of the 
energy of and is concentrated in, respectively, B\ and B2. Outside the balls B\ 
and B2 the asymptotic analysis of the previous section applies. Suppose that the leading 
multipole of is a 2 M -pole and the leading multipole of is a 2 Ar -pole. Denoting the 
radii of B\ and B2 by D\ and D2 we have 

U [1 \x) ~ 1 + u M (x) for x & B x (2.15) 

and 

U {2 \x) ~ 1 + v N (x) for x (£ B 2 , (2.16) 
where um is of the form (|2.9|) and analogously 

v N {x)=ir a ^ oN + l QNn r^+i • ( } 

n=-N 
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Using the translational invariance of the Skyrme energy functional we can assume without 
loss of generality that Bi is centred at X + = (0, 0, R/2) and that B 2 is centred at X_ = 
(0,0, —R/2), where R is so large that B\ and B 2 do not overlap, i.e. R > D\ + D 2 . The 
parameter R will be interpreted as the separation of the Skyrme solitons. Then we define 
the following product configuration 

U R {x) = U {l \x)U {2 \x). (2.18) 

This configuration has degree k + I and finite energy, so that Ur G Ck+i- A lengthy 
calculation performed in J7| shows that 

E[U R ] = Bp®] + E[U®] + AE + O + q (j^—) , (2.19) 

where 

3 

AE = 2J2 d Xl dx 2 (u a M d 3 v a N - v a N d 3 u a M ). (2.20) 

a=l ^3=0 

Now we note that AE is the sum over iso-components of the the scalar interaction terms 
we studied above 

AE a = 2 [ d Xl dx 2 {u a M d 3 v% - v a N d 3 u a M ) = -2V[u a M , v a N }. (2.21) 

Jx 3 =0 

Picking one of the iso-indices, say a = 1, we can use iso-rotations to make sure that the 
first iso-components v} M and v l N are non-vanishing. The result of the previous subsection 
then implies that we can make the multipole interaction energy AE\ non-zero by a spatial 
rotation of the Skyrme soliton with the higher multipole order. 

Now consider the sum 



AE = A Ei + AE 2 + AE 3 . (2.22) 



Following an idea in |]T6| we would like to show that we can always arrange for AE to be 
negative by a suitable iso- rotation of one of the Skyrme solitons. We may assume that, 
possibly after re-labelling the pion fields, 

AE! > AE 2 > AE 3 . (2.23) 

If AE < we are done, so suppose that AE > 0. Since we know that not all AE a vanish 
we can conclude that AE\ > 0. Now perform an iso- rotation of Skyrme soliton 2 by 180 
degrees around the third iso-spin axis. This reverses the sign of v]^ and v 2 N and hence of 
AE\ and AE 2 . The new value of AE is 

AE = -AE 1 - AE 2 + AE 3 

= —AEi - (AE 2 - AE 3 ) < 0, (2.24) 
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since —AEi < and, with our ordering, — (AE 2 — AE3) < 0. 

Thus, the contribution AE to the interaction energy of two Skyrme solitons U^ 1 ' and 
can always be made negative by suitable relative rotations and iso-rotations of the Skyrme 
solitons. Returning to the expression fl2.19|) for the interaction energy and noting from 



fl2~TBl ) that AE falls off like r-W+m+i)^ we conclude that AE is the leading contribution 
for large separation R provided the orders of the leading multipoles do not differ by more 
than 2, i.e. provided that \N — M\ < 2. In that case we can therefore always arrange 
for the interaction energy to be negative for sufficiently large separation. In symbols, for 
sufficiently large R and the appropriate orientations we have 

E[U R ] < E[U W ] + E[U m ]. (2.25) 



3 The existence of Skyrme solitons 

The results of the previous section are not sufficient to derive the existence minima in a 
general sector C k of the Skyrme model. Nonetheless it is instructive to see how far one can 
get with the following two extra assumptions. 

1. When minima exist they satisfy the technical assumptions made in |7j]. In particular 
the field U and the currents Lj are Holder continuous. 

2. The leading multipole of any Skyrme soliton is at most an octupole (order 3). 

The second assumption seems particularly restrictive and deserves a comment. Since 
monopoles (order 0) cannot arise in Skyrme solitons and dipoles (order 1) are known to 
arise in some, the restriction to leading multipoles of order at most 3 means that any two 
Skyrme solitons we consider in the following have multipoles whose orders do not differ by 
more than 2. As mentioned earlier, none of the numerically studied Skyrme soliton violate 
assumption 2. 

Recall from section 1 that Skyrme solitons have rigorously been shown to exist in the 
topological sectors C\ and C-\. Although it is not clear whether the minimum is of the 
hedgehog form ( |1 . 1 1|) , we can use the numerically computed energy of the Skyrmion to give 



an upper bound on the infimum I\. To avoid the discussion of numerical accuracy we use 
the energy of the instanton generated hedgehog field [18| I\ = 1.24 x 127r 2 to estimate 

h < - 12tt 2 . (3.1) 
4 



Using Esteban's weak inequality ( |1.12|) we deduce 
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I k <-12ir 2 \k\ (3.2) 
for all fceZ. Combining this result with the energy bound ( |1.9| ) we also deduce that 

J*<5>. (3-3) 
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for all k G Z and l s G Z, s G N satisfying £) a Z s = ^ an d |Z S | > Thus we can 

weaken the assumption in Esteban's theorem about the existence of minima. The existence 



of a minimum in sector Ck follows if the strict inequality ( 1.13 ) is satisfied for all k G Z and 
alH G Z — {k, 0} satisfying 

1/1 + \k-l\ < -\k\. (3.4) 
4 

Combining this result with our proof of attraction between Skyrme solitons and the as- 
sumptions made above we can prove the existence of minima in the sectors Ck for low values 
of k. We prove the claim by induction and first consider positive k. We know from that 
the minimum exists for k = 1. Suppose that minima exist in C\ with I = 1,2, ...,k — 1. 
Then apply the product ansatz Ur and the inequality (2.25|) to pairs of Skyrme solitons 



U (1) and U {2) of de gree I and k — I for I = 1, ...,k — 1. By assumption the energies of these 
Skyrme solitons are equal to the infima I\ and Ik-i- Then, by definition of the infimum Ik 
AO]) and from (gjg ) 



h < E[U R ] < E[TJM] + E[U^} = k + J (fe _ (3.5) 

for all Z = 1, A; — 1. As long as k < 8 the only integers / satisfying the inequality (|3.4|) are 
in the range I = I, k — 1, so that we have satisfied the condition of Esteban's theorem 
and can conclude that the infimum in Ck is attained. The claim for negative k follows by 
applying the energy preserving reflection map U(x) i— > U'(—x). 

For k > 8 the step in the inductive proof fails because we then need to consider the 
inequality ( |1.13| ) for negative I, too. For example for k — 9, the condition ( |3.4j ) is satisfied 
for / = —1 and k — I — 10. However, since we do not know whether Skyrme solitons exist in 
the sector C 10 we cannot establish the inequality Ig < 7_i+/io with our method. Numerical 
evidence and physical intuition suggests that the inequality Q1.13Q should certainly hold 



when k is positive and I or (k — I) is negative. If it was violated in that range it would 
be energetically favourable for matter (solitons of positive degree) and antimatter (solitons 
of negative degree) to coexist in some sectors of the Skyrme model. However, while this 
seems unlikely, we cannot disprove the possibility with our methods. If one could rule it 
out by other methods, the inductive argument given above could be extended to prove the 
existence of minima in all sectors of the Skyrme model. 

To end, I point out that it seems plausible that our results can also be used to give a more 
direct proof of the existence of minima in all sectors of the Skyrme model. The idea of 
such a proof is that the infimum in one of the sectors, say Cj-, can only fail to be attained 
if it is energetically favourable for configurations in that sector to split into well-separated 
Skyrme solitons. The infimum Ik would then be realised by a virtual solution, made up 
of infinitely separated Skyrme solitons. The number of Skyrme solitons in this virtual 
solution is necessarily finite. For suppose the degrees of the Skyrme solitons are l s , s G N 
with l s = k. By the inequality ( |3.3j ) we can only have the equality 



E4 



(3.6) 
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if S s \ls\ < which means that all but a finite number of the l s must be zero. Suppose 
that l s ^ for s = 1,2. ..S and l s = for s > S. Physically, one expects the following 
dichotomy to hold. 

(a) For given k the infimum 1%. is attained by a Skyrme soliton of degree k. 

(b) The infimum is equal to a finite sum 

s 

= ( 3 - 7 ) 

of energies E\ a = Ii s of S Skyrme solitons of degrees l s G Z — {0, k}. 

With our assumption on the leading multipoles our earlier proof of attractive forces in 
the Skyrme model implies that we can always lower the energy of a finite number of 
infinitely separated Skyrme solitons by chosing appropriate orientations for two of the 
Skyrme solitons and bringing them closer together. Thus the infimum Ik cannot be of the 
form ( |3.7| ) and if the dichotomy is valid, option (a) must hold. 

4 Conclusion 

In this talk I reviewed the result, obtained jointly with Manton and Singer, that arbitrary 
Skyrme solitons attract for suitable relative orientation and iso-orientation, provided the 
orders of their leading multipole moments do not differ by more than two. I also indicated 
how this result, combined with the work of Esteban, can be used to shed light on the 
question of minima in the Skyrme model. However, the approaches sketched here also show 
how difficult it is to establish rigorous results on minima. In particular the assumption on 
the multipoles of Skyrme solitons severely restricts the generality of the sort of argument 
outlined above. 
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